Trap-free manipulation in the Landau-Zener system 
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The analysis of traps, i.e., locally but not globally optimal controls, for quantum control systems 
has attracted a great interest in recent years. The central problem that has been remained open is to 
demonstrate for a given system either existence or absence of traps. We prove the absence of traps 
and hence completely solve this problem for the important tasks of unconstrained manipulation of 
the transition probability and unitary gate generation in the Landau-Zener system — a system with 
a wide range of applications across physics, chemistry and biochemistry. This finding provides the 
first example of a controlled quantum system which is completely free of traps. We also discuss the 
impact of laboratory constraints due to decoherence, noise in the control pulse, and restrictions on 
the available controls which when being sufficiently severe can produce traps. 
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I. INTRODUCTION 

Manipulation by atomic and molecular systems is an 
important branch of modern science with applications 
ranging from optimal laser driven population transfer in 
atomic systems out to laser assisted control of chemical 
reactions [T]. High interest is directed towards control 
of the Landau-Zener (LZ) system — a two-state quantum 
system whose unitary evolution under the action of the 
control e{t) (e.g., shaped laser field) is governed by the 
equation 



ut 



-i(Aa,-|-e(t)(7,)[/f, Ul 



(1) 



where A > 0, cr^: and ct^ are the Pauli matrices. The 
case e(t) = et with constant e was studied by Landau, 
Zener, Stiickelberg, and Majorana [2]. This system has 
been widely applied in physics, chemistry, and biochem- 
istry, e.g., for describing transfer of charge along with its 
energy [3] , photosynthesis [4] , atomic and molecular colli- 
sions, processes in plasma physics fS], Bose-Einstein con- 
densate [7], experimental realizations of qubits, etc. [HllBl" 

Controlled manipulation by a quantum system can be 
formulated as finding global maxima of a suitable objec- 
tive J(e) associated to the system. For example, maxi- 
mizing the probability of transition from the initial state 
|i) to a target final state |f) at a final time T can be 
described by maximizing J{£) — Pj_^£ = |(f|C/^|i)p. A 
control which attains a local maximum of J can be found 
either numerically using the model of the system or ex- 
perimentally. In both circumstances, the first step of a 
common procedure is to apply a trial pulse eg and obtain 
the outcome J(eo), either numerically or measuring it in 
the laboratory. The second step is to make various small 



modifications of eg a-nd find ei which produces maximum 
increase in J . Then ei is used as a new trial pulse and 
the procedure is repeated until no significant increase is 
produced or a maximum number of iterations is reached. 
Of crucial practical importance is to know whether 
J(e) has traps, i.e. local maxima with the values less than 
the global maximum, as necessary to properly choose be- 
tween local (e.g., gradient) and global optimization meth- 
ods [HHIl]- Traps can strongly influence on both theo- 
retical and experimental quantum control studies — they 
determine the level of difficulty of controlling the system 
and can significantly slow down or even completely pre- 
vent finding globally optimal controls. Whereas the anal- 
ysis of traps in manipulation by quantum systems has at- 
tracted high attention [5DH5S] , no examples of trap-free 
quantum systems have been known. Only partial the- 
oretical results have been obtained stating the absence 
of traps at special regular controls. This finding does 
not at all exclude the absence of traps that makes the 
problem open since even a single trap may produce sig- 
nificant difficulties for the optimization if it has a large 
attracting domain j29j . In this work we show that the LZ 
system is trap-free and hence, for example, the only ex- 
trema of J{£) = P^^f for this system are global maxima 
and minima. This finding provides the first example of 
a trap-free quantum control system where unconstrained 
local manipulations are always sufficient to find best con- 
trol pulses. Sufficiently strong constraints on the controls 
may destroy this property and in the end we discuss pos- 
sible limitations for the analysis due to decoherence, noise 
in the control pulses and limited tunability of the control 
strength and time scales in laboratory experiments. 



II. TRAPS AND CONTROL LANDSCAPES 
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Formally, a control field e(t) is a trap for the objec- 
tive J(e) if it is a local maximum, i.e., a maximum with 
the value less than the global maximum, J{£) < Jmax = 



max J(e) (in this work we consider as control goal maxi- 

e 

mizing the objective; if the goal is to minimize the objec- 
tive then traps are local minima). Answering the ques- 
tion whether traps exist for a given control problem is 
crucial for determining proper algorithms and our abili- 
ties for finding optimal control fields. In the absence of 
traps, local search algorithms should generally be able 
to find globally optimal controls (exceptions may occur 
if the initial control is chosen exactly at a saddle point, 
where the gradient of the objective is zero). If the ob- 
jective has traps (perhaps even a single trap with large 
attracting domain) then local search procedures may con- 
verge to local maxima instead of attaining a desired glob- 
ally optimal control and more sophisticated global search 
methods should be exploited for a successful optimiza- 
tion. 

Traps are critical points, i.e., the gradient VJ^ = at 
any trap. Critical points for control objectives J(e) = 
P[^f were studied in seminal works [20', where the ab- 
sence of traps was suggested. The suggestion was drawn 
from the proof that any function of the form f{U) = 
Tt[UpoU^O] (po is a positive matrix and O is Hermitian) 
defined on the unitary group U(n), where n is the system 
dimension, has as extrema only global maxima, global 
minima, and saddles and has no traps. (Extrema of trace 
functions over unitary and orthogonal groups were stud- 
ied in other contexts by J. von Neumann [5IIJ, R. Brock- 
ett [SI], S. Glaser et al [32j, etc.) Then, under the con- 
trollability condition which assumes that any U G U{n) 
can be generated by some control, this result was used to 
conclude the absence of traps for the underlying objective 
functional J{e). Later it was shown that the conclusion 
of the absence of traps requires an additional assumption 
that the map x : e — >■ U^ is non-degenerate |21j , meaning 
that arbitrary infinitesimal variations of e produce vari- 
ations of f/j. in all directions on U{n) [33l[34]. While the 
controllability condition is relatively easy to verify [35], 
checking the non-degeneracy assumption turned out to 
be a hard problem. Moreover, critical controls violating 
this assumption were found [211 El], and even second- 
order traps — critical controls which are not global max- 
ima and where the Hessian H = 5'^J/{de)'^ is negative 
semidefinite were shown to exist under rather general as- 
sumptions |26| . (Second-order traps are not necessarily 
local maxima but effectively they are traps for local algo- 
rithms exploiting at most second order local information 
about the objective; see Chapter 20 of [3B] for a gen- 
eral discussion of the non-degeneracy and second order 
optimality conditions.) These findings led to reconsid- 
eration of the conclusion of absence of traps. Some nu- 
merical simulations suggested that the condition of non- 
degeneracy might be generally satisfied or at least its vi- 
olation does not produce multiple traps |27j, while other 
indicated possible trapping behavior [53J[55]. However, 
numerical search is limited and the extent to which these 
runs span the full space of quantum control possibilities 
is questionable |29| . Hence the problem of proving either 
existence or absence of traps has been remained open. 



III. ABSENCE OF TRAPS FOR THE 
LANDAU-ZENER SYSTEM 

Our main result is that the only critical points of any 
objective of the form J(e) = f{U^), where f/f. satis- 
fies M and f{U) is any function on the special unitary 
group SU(2) which has no local extrema, are global max- 
ima, global minima, and the zero control field e(t) = 0. 
Important examples of such objectives include 

• Transition probability 

This objective is maximized by a control which 
completely transfers the initial state |i) into the de- 
sired final state |f). 

• Expectation of a system observable O 

Jo{e) = Trp^poUpO] 

Here O is a Hermitian matrix representing the ob- 
servable and po is the initial system density matrix. 
The objective is maximized by a control which max- 
imizes quantum-mechanical average of O at time T. 

• Generation of a unitary process W 

Jwie) = l\TriW^U^)\' 

Here W is the unitary matrix representing a desired 
system evolution or a desired quantum gate, for 
example Hadamard gate. Maximum of this objec- 
tive is achieved by a control such that Ut = e^'^W, 
where </> is arbitrary (generally unphysical) phase. 
Factor 1/4 is chosen to have max Jwis) — 1. 

Proof of the main result. We will consider first 
J(e) = J[^f{e). For brevity, we will sometimes omit 
the superscript e in U^ and C/y, and without loss of gen- 
erality set A = 1. Gradient of J(e) = \{i\U^\i)\^ for the 
LZ system has the form |26| 



VJ,(i) = 2^({i\ul\i){i\UTUl<j,Ut\i) 



(2) 



It can be written as VJ^{t) = L{U}azU^) = l{t), where 
L : su(2) — > M is the linear map on the Lie algebra of 
traceless Hermitian 2x2 matrices defined by L{A) = 
23[(i|C/J,|f)(f|[/y^|i)], and l{t) is a real-valued function. 
If £ is a critical control field, then l{t) = and therefore, 
in particular, l'{t) = l"(t) = 0. These derivatives can be 
computed to be 

I'it) - L{-iUJ[a,+eit)a,,a,]Ut) = -2L{u}ayUt) 
l"{t) = -2L{-iUl[<j^+e{t)a,,ay\Ut) 

= -AL{u}a,Ut) + Ae{t)L{u}a,Ut) 



Thus the condition I" - 
e{t) y^ takes the form 



I' = I = ioT any t such that 



L{Uja,Ut) = m^ayUt) = L{Uja,U,) = (3) 

The matrices U^axU^, Ufayllf, U^azU^ are hnearly inde- 
pendent traceless Hermitian 2x2 matrices. They form 
a basis of su(2) and hence (Is]) imphes L{A) = for any 
A G su(2). 

Let |ij^) be the state which is orthogonal to |i). Taking 
A = |i>(i^| + |ii)(i| and A' = *(|i)(i^| - |i^)(i|) gives 

L{A) = => 3((i|4|f)(f|C/j,|ii)) = 
L{A') ^0^ n({i\ul\i){i\Uj,\ij_)) = 

Thus {i\ul\i){i\UT\i±) = 0, i.e. either (i|t/|.|f) = or 
(f |C/T|i±) ~ 0. The former case corresponds to the global 
minimum of the objective (J = 0) and the latter to its 
global maximum (J = 1). These are the only allowed 
critical controls except of e{t) = 0. This finishes the 
proof of the main result for J[_^f{e). 

The analysis above immediately implies that if a lin- 
ear map L : 5u(2) — ^ K satisfies L{UjazUf) = then 
L = 0. Now we will show that it means that the map 
X : e — 7- U^ is non-degenerate everywhere outside of 
e{t) = 0. Since we consider objectives produced by func- 
tions on SU{2) which therefore invariant with respect to 
the overall phase of C/^, we can identify U^ with the cor- 
responding element of SU{2). Small variations around 
[/j. can be represented as C/^ = C/^e*™ « C/j.(l -f Sw), 
where Sw = — i/g UjazUf_Se{t)dt. For the map x to 
be non-degenerate, [/^ should span a neighborhood of 
U^ that in turn requires Sw to span su(2). If Sw does 
not span su(2), then there exists A £ su(2),yl ^ 
such that {A,Sw) = Tr(^^(5w) — for all Se and hence 
LAiUlazUt) := TiiA^fUJazUt) = 0. This is possible only 
if A = and hence the map can not be degenerate. 
Therefore our result immediately implies the absence of 
traps at any e ^ for any objective functional J{Uj^) 
which has no traps if considered as a function on SU{2). 
This includes important objectives Jo = Tr[U^pQU^O] 
for maximizing expectation of a system observable O 
and Jw — (l/4)|Tr(Vl^1'C/j.)p for optimal generation 
of a unitary process W (e.g., for unitary gate genera- 
tion). These objectives appear to be trap-free for the 
LZ system since functions fo{U) = Tt\UpoU'^0] and 
fw{U) = {1/4:)\Tt{W^U)\^ have no local maxima on 
SU(2) ^. 

The control e{t) = requires a separate consideration 
since the condition l"{t) = for e = can not be used to 
conclude L{Uj axU^) — 0. This control is however not a 
trap for example for Jj_j.f as shown by direct computation 
in the Appendix. 
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FIG. 1: (Color online) The control landscape of Jo^i(ai, 02) 
for the LZ system controlled by piecewise constant controls 
(A'' = 2, T — 10, A = 1). The landscape possesses multiple 
traps (local maxima). 



IV. DISCUSSION 

Now we discuss important limitations for the present 
analysis. No real-world system will perfectly evolve ac- 
cording to Eq. (fTl) and three general kinds of deviations 
from the ideal situation include decoherence effects, de- 
viations of the actual control from the intended one due 
to noise or imperfections of the laboratory setup, and 
limited tunability of the control strength and time scales 
in laboratory experiments. While we consider the sys- 
tem as evolving according to the Schrodinger equation 
with unitary evolution, in real circumstances it can ex- 
perience additional influence of the environment which 
causes the dynamics to be non-unitary. We also assume 
that any shape of the control e{t) is available, whereas 
typical pulses are either piecewise constant or finite sums 
of cosines and sines at certain fixed frequencies. These 
assumptions are common for the first step of control land- 
scape analysis which deals with the ideal situation of 
noiseless unconstrained controls. The next step upon es- 
tablishment of the ideal landscape properties is to study 
the effects of possible deviations, which we discuss below 
for the LZ system. 

The requirements on the available control fields (e.g. 
on their strength and time scale) necessary for the con- 
clusion of the absence of traps for attaining maximal ob- 
jective value are such that the available controls are suffi- 
cient to guarantee controllability of the system. Minimal 
control time for the LZ system can be estimated using the 
fundamental theory of optimal control at the quantum 
speed hmit as Tqsl ~ Ai^p"^ arccos(|(i, f)|), where A_Eo 
is the energy variance of the free Hamiltonian Hq = A-a^ 
calculated on the initial state [TDj. Hence our analysis 
applies to any final time T > ttAEq^. If for a given 
physical system decoherence effects occur on a time scale 
slower than AEq^ , they can be neglected when the con- 
trol is implemented in the time optimal fashion. This 
shows that while finite-time [37] and decoherence [55HiD] 
effects can be important for the LZ system, they do not 
modify the trap-free landscape property as soon as final 
time T is sufficiently smaller that the relaxation time and 
at the same time is not too small to violate controllability 




FIG. 2: (Color online) Probability of trapping as a func- 
tion of TV for piecewise constant controls (T = 10, A — 1). 
For every point, 10'^ runs of MATLAB realization of the 
Broyden-Fletcher-Goldfarb-Shanno (BFGS) optimization al- 
gorithm where performed each starting at a random initial 
control a — (ai, . . . ,ajv) [41]. Initial control amplitudes are 
uniformly distributed in the range Ui £ [A, A] but are allowed 
to escape this range during the search. The search is defined 
as trapped if the attained objective is less than 0.99. Trap- 
ping may occur due to the presence of local maxima and/or 
principal impossibility of attaining the objective value greater 
that 0.99 with available controls. Probability of trapping is 
estimated as a fraction of trapped runs among all 10'' runs. 



of the systems. 

An extensive numerical analysis of control landscapes 
for multi-level model systems with realistic laboratory 
control fields is provided in [57]. To analyze the role of 
limitations on the available control fields for the LZ sys- 
tem, we numerically estimate the probability of trapping 
when available controls are piecewise constant controls of 

the form e(t) = Ej=i aiX[t,, ti+i](i), where X[ti,ti+i](0 = 1 
if t G [ii,^i+i] and zero otherwise and a^ are the control 
parameters. Typically, N « 100 and control amplitudes 
are constrained within certain ranges, say Oj e [—A, A]. 
Exact solution for piecewise constant controls can be ob- 
tained for example in the simplest case A^ = 1. The 
objective for maximizing the probability of spin flip 
Jq^i — |(0|[/^|l)p by a constant control e{t) = a can 
be computed to be Jo^iia) = sin^(TVl + a'^)/{l + a^). 
Its traps (local maxima) are given by solutions of the 
equation tan(Tv'l + a'^) = T^/l + a^; the correspond- 
ing objective values are Jo->i(a) = T^/{1 + T^ + T'^a). 
Control landscapes for a two-dimensional control space 
{N = 2) are more complex. Fig. [I] shows as an exam- 
ple the control landscape of Jo->.i(ai, 0,2)- The landscape 
has multiple local maxima showing that significant re- 
strictions on the control space in an originally trap-free 
system may produce traps. Fig. [5] provides the numeri- 
cally estimated probability of trapping for piecewise con- 
stant controls as a function of N. The probability of 
trapping becomes negligible already for N — 10-15 that 
means that limitations on the number of components N 
of available laboratory control fields have minor effect 
already for iV > 10 and hence should be negligible for 



realistic case N w 100. 

In the laboratory, actual controls may deviate from the 
designed numerically optimal pulse due to noise and im- 
perfections of experimental setup [35J . These noise effects 
can influence on the landscape structure by decreasing 
the maximal objective value. We adopt the general the- 
ory of [13] to analyze this influence for the LZ system. 
Let Soit) be an optimal control in the ideal situation of 
absence of noise. In the presence of a random noise (,{t), 
the actual control will fluctuate as e{t) = eo{t) + g{t)£,{t), 
where g{t) = 1 for additive noise and g(t) ~ £o(0 for 
multiplicative noise. A weak noise modifies the averaged 
objective as 



T T 



E[J(eo)] « JieoY 



where ^{^{t.t') 



Ho{t,t')g{t)Q{t')¥.[i{t)S,{t')]dtdt' 







5'J 



is the Hessian of the ob- 



5eo(t)6eo{t' 

jective computed at the optimal control field Eq ^^nd 
E[^(i)^(i')] is the autocorrelation function of the noise. 
Since Hessian is negative semidefinite at the maximum, 
the noise generally decreases the average fidelity. The 
objective for additive (AWN) and multiplicative (MWN) 
white noise with autocorrelation function E[^(i)^(i')] = 
aS{t — t'), where a^ is the variance of the noise amplitude 
distribution, takes the forms 



EAWN[J(eo)] « Ji£Q) + Y I ^oit,t)dt 



Emwn[J(£o)] ~ J{£o) + Y I 'Ho{t,t)\so{t)\'dt 



The last term in these equations is the noise-induced de- 
crease —'D{eQ,(T,T) of the objective (such that E(J) « 
J{eo) - V with X> > 0). The diagonal of the Hes- 



sian for J 



Ji_ 



/j_^f can be shown to be 'Ho{t,t) = 

Therefore 
I?(eo, o", T) for J = J[^f is majorized by 



rt. 



'2\{i\Ula,Ut\iA_)\^ so that \'Ho{t,t)\ < 2. 



2?AWN(eo,o-,r) < (T^T 
2?MWN(eo,o-,r) < (T^E 



where E = L \eo{t)\'^dt is the total energy of the pulse. 
The diagonal of the Hessian for the objective Jw is 
HQ{t,t) = —2 and therefore for this objective Pawn — 
a^T and I'mwn = (J^E. It then follows that in both cases 
the influence of a weak AWN can be minimized by us- 
ing time optimal controls, while minimizing weak MWN 
can be done by selecting less energetic pulses among all 
optimal pulses. 

If the ideal landscape has multiple global optima with 
different Ho(t, t), then the noise induced decrease of ob- 
jective can be different at different optima that can pro- 
duce traps in the non-ideal landscape. Weak decoherence 



operates similarly to weak noise and can also produce 
traps in the ideally trap- free landscape [20]. These de- 
viations from the ideal situation should be avoided to 
reveal the trap-free landscape property by either oper- 
ating in the time optimal regime or using weak optimal 
controls to combat MWN. Strong noise and strong deco- 
herence that can significantly modify the landscape are 
out of scope of this discussion. 

Conclusions. — This work shows that unconstrained 
manipulation in the Landau-Zener system is free of traps 
and hence unconstrained local search for optimal con- 
trols is always able to find best optima. The impact on 
this result of laboratory limitations due to decoherence, 
noise in the actual control pulses, and restrictions on the 
available control fields is discussed. 



in Se as 

Wt = I + Ai + A2 + oi\\6e\\^) 

Ai = -i I dt5e{t)Vu 
Jo 

A2 = ~ [ dh r dt2Se{h)5eit2)Vt,Vt, 
Jo _ JO 
This gives the perturbation expansion for the objective 

(here |f') =ei^'^-|f)) 

J{Se) = \{i'\I + Ai+A2 + ...\i)\' 

= |(f'|i)|2 + JJi(fe) + SJ2iSe) + o(||fc||2) 

where 
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Appendix 



SJiiSe) = 23?((f'|i)(f'|Ai|i)) 

SJ2{Se) = |(f'|Ai|i)|2 + 23?((f^(f'|A2|i» 

Hence the variation of the objective satisfies (note that 
|(f'|i)P = J(0)) 

6 J = J(fc) - J(0) = 6Ji{Se) + 5J2{5e) -f o{\\5ef) 

If e = is a critical control, then 5Ji [Se) = for any 
Se. We will show the existence of controls Sei and 5e2 
such that 5J2{5ei) and 5J2[Se2) have opposite signs. It 
is sufficient to choose fci and 5£2 to satisfy (f'|Ai|i) — 0, 
e.g. 



T ^T 

dtde,{t) cos2t= dtSei{t) iim2t ^ 0, i = l,2. 



Here we prove that the control e{t) = is not a trap for 
state-to-state transfer described by the objective J{e) = 
Jj_^f(e). The evolution operator produced by e{t) — 

has the form Ut = e"**'^-. Therefore Vt := UJa^U^ ^ Since Vt^Vt^ = cos2(ii - is) + icr^ sin2(ti - tj), we have 
cos(2i)cr2 -I- sin(2i)cry and the gradient of the objective is 



VJe=oit) = cos2t ■ L{a^) + sm2t ■ L{ay) 

If e(i) = is a critical point, then WJe=o{t) — for any 
t € [0,T], and hence L{az) = L{ay) = 0. If a := L(o-^) = 
0, then L = on su(2) and similarly to the proof of the 
main result we conclude that e = is not a trap. 

Now consider the case a ^ 0. In this case |i) and |f) 
are such that e = is neither a global maximum nor 
global minimum. The evolution operator produced by a 
small variation of the control 5e can be represented as 



U^^ 



e '-"^^Wy, where Wt satisfies 



Wt ^ -i5e{t)VtWu Wo=I 
The operator Wt can be computed up to the second order 



5J2[5e) = 2/" dh I dt2Se{ti)Se{t2)(j{0)cos2{ti-t2) 



+asm2{ti — i2) 

2a dti I dt2Seiti)6e{t2)sm2{ti-t2) 



Assuming for simplicity that T > tt, we take fci(t) — 
Xlo,^]{t) and fe2(i) = cos{'it)x[o,Tr]{t), where X[o,7r]it) is 
the characteristic function of the interval [0,7r]. Then 
SJ2{Sei) = —TTa and SJ2{Se2) = ira/d. Therefore for 
a ^ there exist control variations around e{t) = 
increasing the objective and control variations decreasing 
it. This implies that e(t) = is neither a local maximum 
nor minimum. 
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